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Abstract 

The liquid-gas phase transition in strange hadronic matter is studied utihzing an extended 
Furnstahl-Serot-Tang model with nucleons and hyperons. The system is treated as of two compo- 
nents. The phase transition is analyzed by investigating the stability of the system and Gibbs con- 
ditions for phase equilibrium. A two-dimensional binodal surface resulted from the two-phase equi- 
librium, namely the phase boundary, is obtained. For each temperature ranging from T = 8MeV 
to T = 12MeV, a limit pressure on the binodal surface section is found, while a critical point is 
spotted for the temperature around T = 13MeV. The Maxwell constructions are also illustrated 
to give a vivid description of the course of the phase transition. Moreover, the entropy per baryon 
and heat capacity per baryon as functions of temperature are examined. The entropy is continu- 
ous during the phase transition but the heat capacity is discontinuous, indicating that the phase 
transition is of second order. By these efforts, the L-G phase transition can be concluded to exist 
in the strange hadronic matter. 
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I. INTRODUCTION 



Since the nuclear interaction is very similar to the Van der Waals potential between 
molecules, the possible existence of liquid-ffas(L-G) phase transition in nuclear matter 
was predicated by theoretical studies l|, S Hi • This predication was supported later by 
intermediate-energy heavy- ion collisions jj| and high-energy proton-induced reactions 5|. 
The idea of an L-G phase transition in nuclear matter emerged and attracted much 



interest 
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12l |. This interest has increased in the last decade, with the 
attempts by the EOS Collaboration to extract exponents of 1 GeV/nucleon Au nuclei with 
a carbon target jisl and with the extraction by the ALADIN/LAND Collaboration of a 
caloric curve resulting from the fragmentation of the quasiprojectile formed in the colli- 
sion Au-|-Au at 600 MeV/nucleon exhibiting a behavior expected for a first-order L-G phase 
transitional . Very rich information has already been extracted from exclusive experimental 
data during the intermediate energy heavy ion collisions, as the new generation An detectors 
have been developed and are now operating at different accelerator facilities, such as Dubna, 
GANIL, GSI, LNL, LNS, MSU and Texas A-m|i3- 

This L-G phase transition in nuclear matter occurs at lower density than the normal 
nuclear density po 16| . The two phases are of cold nuclear Fermi liquid on the one hand and 
a nuclear gas consisting of free nucleons on the other hand. For symmetric nuclear matter, 
having equal number of protons and neutrons, this L-G phase transition is known of first 
order. However, Miiller and Serot have indicated that the phase transition in asymmetric 
nuclear matter is of second order by Ehrenfest's definition |l7l|. The entropy and the volume 
as the first order derivatives of the chemical potentials are continuous and the heat capacity 
as the second order derivative is discontinuous at the point of phase transition. The phase 



coexistence surface of this binary system of two phases is two dimensional 
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24|. 



Since the first hypernucleus was detected in emulsion by Danysz and Pniewski |25l| in 1953, 
strangeness carried by s-quark has opened a new dimension for the studies in nuclear physics. 
With the hyperons included in, nuclear matter is extended to strange hadronic matter (SHM). 
Employing different models, the equation of state and the stability of the SHM have been 



investigated 
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However, besides the deconfinement phase transition to 



the quark gluon plasma(QGP), whether there exists an L-G phase transition in the SHM 
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is still an open question in theoretical physics. As the interactions between hyperon and 
nucleon or between hyperons are also of the behavior of the Van der Waals potential, it is 
reasonable to guess that there will also be L-G phase transition in the SHM. It is the main 
purpose to study this problem. 

It is believed that the descriptions of nuclear matter and finite nuclei are ultimately 
governed by the physics of low-energy quantum chromodynamics(QCD). Because of the 
nonperturbative properties of QCD, it is very difficult to establish a fundamental theory of 
nuclear system from QCD. One usually adopts various phenomenological models based either 
on hadron degree of freedom or quark degree of freedom. In our previous papersj^, we 
suggested an effective model, constructed by introducing hyperons in the Furnstahl-Serot- 
Tang(FST) modeljs^, Is^, to study the saturation properties and stabilities of the SHM. 



The details of this model can be found in ref. 
In refs 



Here we only give a short description. 



321, we extended the original FST model to include not only nucleons and a. 



uo mesons, but also A, S hyperons. A new hyperon-hyperon interaction mediated by two 
additional strange mesons a* and (j) , which couple only to hyperons, is introduced. The 
reactions A + A— A + A— + and their reverses are considered. In this model, 
we did not consider the mixture of the S hyperons. The reason, as was explained in ref. 31| , 
is twofold. First, the S potential in nuclear matter at the saturation density is rather 



35| to f/s = -25 ± hMeV 



uncertainly predicted, ranging from completely unbound 

A. poised Ba.e, e:alH 

, systems involving E's together with nucleons or A's 

generally will be unstable with respect to the strong decays E+A^ A+A^ or S+A S+A^. 
Secondly, the Q values for the strong transitions H + N —>■ A + A^, S + S A + A, 
S + A->S + A^andS + E:->A + Hare about 78, 156, 50 and SOMeV, respectively 3]- 
To Pauli-block these processes, we need a rather high density of A. On the other hand, 
the Q value ofS + A^— i>A + Ais only about 28MeV. This process can be Pauli-blocked 
by a relatively low A density. Based on this model, the stability, the density, temperature 
and strangeness fraction dependence of the effective masses of nucleons and baryons, the 
pressure, the free energy and the equation of state are studied, and the results are reasonable. 

The paper evolves from an attempt to study the L-G phase transition of SHM by em- 
ploying the extended FST model. In Sec. II the extended FST model is laid out. The 
SHM is modelled into a two-component system ready for the discussion of phase transition. 
The thermodynamic quantities for this model are briefly derived out. By analyzing the 
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equilibrium conditions the phase sections are addressed in Sec. III. The thermodynamical 
properties and the phase transition are discussed and the order of the phase transition is 
also examined. Sec. IV finally summarizes the main results about the existence of the L-G 
phase transition in the SHM. 



II. THE EXTENDED EST MODEL 



The original FST model is extended by including A and S hyperons and two additional 
strange mesons a* and to describe the hyperons and the interactions between them as 
suggested in ref. [31j. Since we are studying the unpolarized system, the vr meson has no 
infiuence on the system and is omitted. To reduce the degrees of freedom, we need some 
assumptions. First, we assume S~ and S° appear in equal amount. Also we have protons 
and neutrons in equal amount, which is similar to protons and neutrons in symmetric nuclear 
matter. These assumptions mean that we are looking at matter with zero isospin. But the 
remarkable new degree of freedom, strangeness is included. The lagrangian of the extended 
FST model is presented as follows: 
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In mean-field approximation, the lagrangian can be written as 

C-MFT = + gujNl^Vo - Mn + gsNcro)'ipN 

+^pA{iYd^ - gu;Al^Vo - g<pAl°<Po - Ma + gsAc^o + fi'a*AO-o)^A 











r 










~ 4. 



(1) 



(2) 



where the meson field operators have been replaced by their mean field values: 0o, Vo , co 
and (Jq. Qij are the coupling constants of baryon j to meson i field. By using the standard 
technique of statistical mechanics, the thermodynamic potential ft is obtained 

^ = m[(l-|)^(^ln(l-|)4).l] 

-2kBT{Y,\^ [1 + e-^(^*(*^)-'^^] + ^ In [1 + e-''(^n'=)+-i)]} (3) 
where (5 is the inverse temperature and V is the volume of the system. 



E*{k) = ^Mf + P (4) 
with the efi^ective masses of the hyperons and nucleons 

M* = Mi- g,iao - g^nc^o (i = A, S) , (5) 

M* = Mi- Qs^ao (^ = N) (6) 

The mean-field values (f)Q,VQ , (Tq and CTq are determined by the corresponding extreme condi- 
tions of the thermodynamic potential. For instance,the equation for meson is determined 

by 

'Q^\cro,Vo,al,iJ.i — (7) 

in explicit form 

"^|0o - g4>KpBA - 9<psPBS = 0. (8) 
The baryon densities pBi is given by 

PBi = {^tA) = ^J dke [m {k) - Ui (k)] (9) 

where g^j = 4 for i = A?" or S, ^fj = 2 for i = A. The baryon and anti-baryon distributions 
are, respectively, expressed as 

mik) = {exp[(3{E*{k) - Ui)] + 1}-' (10) 

and 

n,{k) = {exp[P{E*{k) + u,)] + 1}"^ (11) 
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in which Vi are related to chemical potential /ij by 



IXn — I'm + QunVq, 
Pa = i^A + QljaVo + y4>AVo, 
Pe = i^s + gu^sVo + g^3(t>o- 



(12) 



In the system with equal number of protons and neutrons as well as equal number of and 

the chemical equilibrium condition for the reactions A + A ^ n + S° and A + A ^ p + 
reads 

'^Pa = Pn + Ps- (13) 

Eq.(13) implies that only two components are independent among A^, A and S. The system 
of SHM described by the extended FST model is a two-component system. The nucleons {p 
and n) and the hyperons (A, S~ and H*^) play the role of different components, respectively. 
The strangeness fraction is introduced as 

Pba + 2pBH 



fs 



PB 



(14) 



mm. 



m 



which plays the similar role as that of the asymmetric parameter a = {pn — 
the asymmetric nuclear matter. We can use the same method as that in ref.(l7l.ll 
to address the L-G phase transition. 

Following the usual procedure of statistical physics, we can easily calculate the other 
thermodynamic quantities from thermodynamic potential Q. For example, the pressure and 
entropy density are calculated by formulas p = —Q/V and S/V = —{dVL/d{l/l3))v,fijV = 
{dp/d{l/(3)) 

y,^i,o-o,Vb,0o,o-S- results are expressed as follows 



p 



E 



dk 



Em 
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ni{k) + ni{k)\ - Hq 
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(15) 



9i 



dk- 



■ (3\E*{k)-Ui)eM(^m{k)-Ui))] 



^ J ^■^Et{k) \ [exp {Pm{k) - + 1]^ 
^\El{k) + v,)eMPm{k) + v.;))] \ 



[exp [/5((E*(A;) + z/,))] + ir 



(16) 
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III. L-G PHASE TRANSITION 



In this section, we employ the extended FST model to investigate the L-G phase transition 
in the SHM at different strangeness fractions. In the calculation, the parameter set Tl given 
in ref.jsil is used. 

As the hot and dilute SHM is only obtained on earth by the relativistic heavy ion collision, 
the investigations of the thermodynamic properties of the SHM are based on two basic 
assumptions. One is that we can apply equilibrium thermodynamics for such a small system 
of only few hundred constituents at the most. The other is that a thermalized uniform system 
is formed in heavy ion collision before the multi-fragmentation takes place js^. Although the 
equilibrium analysis oversimplifies the study of the SHM, we still follow the thermodynamic 
approach in the reason that this can give some concrete descriptions of the phase structure 
of the SHM and characterize certain aspects of the evolution. 

We would like to discuss the stability of SHM first. To address the stability, we consider 
the Helmholtz free energy density JF(T, pi) defined by temperature T and baryon densities 
Pi- A system is stable against separation into two phases if its free energy is lower than the 
free energy in all two-phase configurations. This statement is formulated as 

^ (T, p,) < (1 - \)T (T, + (T, pD (17) 

with 

p, = (1 - \)p[ + \pl (18) 

The two phases are denoted by a prime and a double prime. In asymmetric nuclear matter, 
we may choose pi as nuclear density p and asymmetry parameter a. The parameter A = 
V" /V specifies the volume fraction of the phase with double prime. Eq.(18) ensures the 
overall baryon densities are conserved. In the SHM here, we may use variables pB, fs instead 
of Pi. The above stability condition implies the following set of inequalities[17] 



Of'BJrj. \Of'B/TJ, 



P =pM^ >0 (19) 



where p, Pn and p= are, respectively, the pressure, chemical potentials for nucleon and S. 
The first inequality Eq.(19) means that the isothermal compressibility is positive, i.e. the 
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system is mechanically stable. The second condition Eq.(20) reflects the special character of 
the binary system. It expresses " diffusive stability" which guarantees that energy is required 
to change the concentration in a stable system, while holding the remaining variables {p and 
T) fixed. 

28| | and show the possibility for an L- 



Now we begin to discuss the region < 1.0 
G phase transition in the SHM. To be more specific, Fig.l shows the pressure-baryon 
density isotherms at temperature T = 8MeV with different strangeness fractions {fs = 
0.1,0.3,0.5,0.8 and 1.0). One can see that there is always a mechanical unstable section 
where condition (19) is violated for any strangeness fraction fs considered. This means 
that the system always encounters an unstable region and has to separate into two phases 
to maintain the mechanical stability. This phenomenon remains up to T ^ ISMeV. At 
temperature T = ISMeV, as is illustrated in Fig.2, there is a mechanical inflection point 
which satisfles the condition 

dp d^p 



(21) 



dp (9p2 1 fs-fs 

with = 0.5 and the corresponding pressure p = 0.24MeV"/m~^, pB = 0.063/m~^. 
For fs > 0.8, it starts to decrease in high density region and mechanical unstable section 
will occur. However, one cannot flnd the solution which satisfles the Gibbs conditions 
in this region (See below in Fig. 7). To flnd the temperature dependence of f^'^, we also 
show the situation for T = lAMeV in Fig. 3. The mechanical inflection point is around 
p = O.lSMeVfm-^ and pB = OMfm-^ on the isotherm with fs = 0.1 for T = UMeV. We 
flnd the f^^ decreases when the temperature increases. In Fig. 4, we present the pressure- 
density isotherms of flxed strangeness fraction fg = 0.5 at different temperatures. One can 
see again an inflection point on the curve with T = 13MeV. 

Next, we will discuss the chemical instability by showing the chemical potential isobars 
for nucleons and S against strangeness fractions fs at temperature T = 13MeV for pressure 
p = 0.10, 0.18, 0.28 and 0.35MeV/m~^, respectively, in Fig. 5. There is an inflection point on 
the curve with pressure p*^ = 0.28Me\^/m"^. This pressure is called critical pressure 17. 2ol|. 



above which the system is stable but below which condition (20) is violated and the system 
becomes chemical unstable. The critical pressure p'" is determined by the inflection point 
condition: 

^dpN\ fd'pN\ ^22^ 



^fs J T,p \ ^fs J T,p 
8 



The result (jf , f^) = {0.28MeV fm ^, 0.65) defines a critical point for a given temperature 
T = 13MeV. 

When the system separates into two phases, the two coexistent phases are governed by 
the Gibbs conditions, which say that for two phases (or more) the system should be in 
chemical, thermal and mechanical equilibrium. 

/ij (t, p^ f^) = /.^^ (t, p^, /f ) , (g = iV, S), (23) 

/(T,p^/^)=p^(T,p^,/f), (24) 

where the superscripts L, G denote the liquid and gas phases, respectively. Here we only 
consider the situation of two-phase case. As to more than two phases, there is some dis- 
cussion in ref. jl^l- The solution {p^, f^; p^, f^} specifying the two separate phases in equi- 
librium can be easily found through a geometric approach. As an example, we present in 
Fig. 6 the chemical-strangeness fraction isotherms at temperature T = lOMeV and pressure 
p = 0.09MeVfm~^. The desired solutions form the edges of a rectangle shown by dotted 
lines. The phase with lower(higher) fs corresponding to a higher(lower) density is liquid(gas) 
phase. It is obvious that the strangeness fractions are different in the two phases. The col- 
lection of all such pairs {p,T, f^) and {p,T, f^) forms a binodal surface of two dimensions, 
which defines the phase separation boundary. The shape of the whole binodal surface is 
similar to that in Ref. [itI I20I 1. To explain clearly we show in Fig. 7 a section of the binodal 
surface at T = 13MeV. In Fig. 7 the binodal curve is divided into two branches by zero 
strangeness point and critical pressure(CP) point. The branch for larger fs represents the 
gas phase with lower density and the other for the liquid phase with higher density. In an 
isothermal compression, a system with zero strangeness evolves until it encounters the bin- 
odal and then remains there until the transition is completed. In an isothermal compression 
with fs > 0.8, the L-G phase transition cannot take place because the isothermal line does 
not cross the binodal surface. At CP point, the liquid branch and gas branch joint together 
smoothly and then the two phases can no longer be distinguished by their densities. The 
CP point determines the maximum pressure in two-phase region. If the system is repre- 
sented by the point in the section lower than the gas phase branch, the system is in the gas 
phase, whereas if the system is at the point higher than the liquid phase branch, it is in the 
liquid phase. Between the two branches, the system is in the mixed phase with a special 
proportion between liquid and gas phases. The feature of the binodal curve at temperature 



17, But for the 



T = 13MeV is very similar to those in asymmetric nuclear matter 
temperature lower than T = ISMeV, the critical points could not be attained. Instead, we 
can always find a limit pressure in the binodal surface. As an example, we show in Fig. 8 
the section of the binodal surface at the temperature T = lOMeV, where a limit pressure 
Piim = 0.095MeVfm~^ is denoted by a horizontal dotted line. When pressure is higher than 
the limit pressure pum, the L-G phase transition cannot take place. The cut-off behavior 
for the binodal surface has also been found in asymmetric nuclear matter described by FST 
model with density-dependent NNp coupling S'pQ, H]- 9p constant, the cut-off be- 
havior will not happen so that a critical point could always be found and no limit pressure 
is obtained. In this case, the L-G phase transition becomes more complicated. The limit 
pressure exists when T < 13MeV even all couplings are constant. 

The Maxwell construction of the phase transition can narrate a phase transition clearly 
and give a concrete proof for the existence of the L-G phase transition. For a given tempera- 
ture, we will expatiate the behavior of SHM under isothermal compression. At T = lOMeV, 
for instance. Fig. 8 illustrates the situation of the isothermal compression. Assume that the 
system is initially in the gas phase with strangeness fraction fs = 0.15. During the com- 
pression, the phase boundary is encountered at point A. At this point, the liquid phase is 
about to emerge at point B, where the strangeness fraction = 0.05. During the whole 
compression, the total strangeness fraction of the system conserves but the strangeness frac- 
tions for the liquid and gas phases are changeable. As the compression proceeds, the gas 
phase evolves along the gas phase branch of the binodal surface from point A to D, while 
the liquid phase evolves from B to C along the liquid branch. When the pressure of the 
system continues to increase, the system leaves the two-phase region at point C, which has 
the same strangeness fraction with point A. Correspondingly, at point D on the gas branch 
with strangeness fraction fs = 0.4, the gas phase disappears. Since the two points A and 
C are at different branches of the binodal surface, the SHM has undergone an L-G phase 
transition. 

To configure out the evolution of the system between point A and C, we must solve the 
following conservation equations 

PB = {1- \)P% + (25) 

PBfs = {l-\)p%f? + \p'hfs. (26) 
10 



for the baryon density, strangeness fractions on the binodal surface and the proportion A for 
given total baryon density pb and strangeness fraction fg. The proportion A is at point A 
and 1 at point C, and it runs through the interval [0, 1]. The result is the generalized Maxwell 
construction in the binary system. The isotherms corresponding to Fig. 8 are drawn in Fig.9. 
The system dose not evolve along the unphysical curve(the dotted line). It follows the nearly 
strait solid line between point A and C, which is the result of the Maxwell construction. Each 
point on the line corresponds to the stable configuration at each intermediate density during 
the phase transition. Moreover, if the system is initially in the gas phase with fg — 0.5, the 
process of isothermal compression will make that the system begins at the gas phase, then 
enters a two phase region and becomes unstable at the limit pressure pum- 

Finally, wc would like to discuss the order of the L-G phase transition in the SHM. By 
Ehrenfest's definition, the first order phase transition is characterized by the discontinuities 
of the first order derivatives of the chemical potential, such as the discontinuities of entropy 
and volume, while the second order phase transition unfolds the discontinuous behavior for 
the second order derivatives of the chemical potential, such as the heat capacity. Using the 
entropy density of the SHM in Eq.(16), we have the entropy per baryon as 

.(T,p,/,) = ^^i^^. (27) 
Pb 

The total entropy per baryon of the system is calculated by the equation: 

s = (1 - A)s^ + Xs^. (28) 

The entropy per baryon as a function of temperature for a matter of fg = 0.50 during the 
phase transition at a fixed pressure p — 0.17MeV fm~^ is delineated in Fig. 10. It is obvious 
that the entropy per baryon varies continuously during the phase transition. To make our 
result more transparent, we calculate the heat capacity Cp of the SHM. 



TV^ [dk ^ 
4" 67r2 J 



4 



E*{k) 



+ 



dT'^ dT^ 



(29) 



where 



d^rii 



j{E*{k)-Ui) 



{E*{k) - Vif{e'^i^*i^^^-''') - 1) - 2T{E*{k) - yi){e^i^i^''^-''') + i) 



(30) 
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{E*{k) + t/i)2(e^(^*(*^)+"0 - 1) - 2T{E*{k) + 1/0(6^^'''^'=^+"') + 1) 



(31) 



c 

The heat capacity per baryon is expressed as Cp = In Fig. 11. the specific heat capacity as 
a function of temperature for a fixed pressure p = 0.17MeV fm~^ for fs = 0.50 is presented. 
A finite discontinuity of the heat capacity is observed, which demonstrates clearly that the 
L-G phase transition in SHM is of second order. 

Till now, we can draw out a definite conclusion that an L-G phase transition do exist in 
the SHM and it is of second order. 



IV. SUMMARY 



In this paper we have employed an extended FST model with nucleons and A, S hyper- 
ons to describe the thermodynamical properties of SHM. By using the reactions between 
hyperons and nucleons and making some assumptions, we have simplified the SHM into a 
two-component system, which is mathematically feasible to be discussed. When having set 
the strangeness fraction to be zero, the model can also obtain results, which match that of 
symmetric nuclear matter. Applying the model, we have investigated the L-G phase transi- 
tion of the SHM. The mechanical unstable region always exists for all strangeness fractions 
at the temperatures lower than T — ISMeV. Meanwhile, chemical instability is also found 
in chemical potential-strangeness fraction isobars at these temperatures. By using the Gibbs 
conditions of phase equilibrium, we have figured out the binodal surface as the phase sep- 
aration boundary and have found a limit pressure at each temperature, above which there 
is no L-G phase transition. The limit pressure gives rise to a cut-off on the section of the 
binodal surface. But for temperature as high as T = ISMeV, a critical pressure has been 
obtained in the binodal surface section. At the critical pressure point the liquid and gas 
branches joint together, resulting in a closed binodal curve. The Maxwell constructions for 
the phase transition have been depicted in detail. The entropy and heat capacity have also 
been examined and the L-G phase transition in the SHM has been determined to be of 
second order. 
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FIG. 1: Pressure as a function of baryon density at temperature T = SMeF for various strangeness 
fractions fg- 



15 




16 




17 



0.6- 



0.4- 



-0.4 



1 


1 .. . 1 / /I / /I 
14/ / / / 






/l3 / / / / 






/ / ^ ^ / / / 






/ / / ^ y / / 






y / / / 






— ^/ 1 






^^^^ 8 / 




f =0.5 

s 







0.00 0.05 0.10 0.15 0.20 0.25 0.30 



Fig.4 



FIG. 4: Pressure as a function of baryon density for fixed strangeness fraction fs = 0.5 at various 
temperature. A mechanical inflection point is at isotherm T = ISMeV, below which, isotherms 
have the mechanical unstable regions. 
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FIG. 5: Chemical potential isobars as a function of the strangeness fraction fs at fixed temperature 
T = 13MeV. The curves labeled a through d have pressures p = 0.10, 0.18, 0.28 and O.S5MeVfm~^ 
respectively. Point A denotes the inflection point at the critical pressure p'-^ = 0.28MeV fm~^. 
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FIG. 6: Geometrical construction used to obtain the strangeness fraction and chemical potentials 
in the two coexisting phases at fixed temperature T = lOMeV and pressure p = O.OQMeV fm~^. 
The two solid crossing curves for iin (or are the different parts of the same continuous curve. 
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FIG. 7: Binodal surface at T = ISMeV. A critical pressure(CP) and the zero strangeness are 
indicated. The upper branch is for hquid phase and the lower for gas phase. Two branches joint 
at CP point. 
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FIG. 8: Binodal surface at T = lOMeV. The binodal surface is cut off at a limit pressure 
p = 0.095MeVfm~^. The point A through D denote phases participating in a normal phase 
transition. Another phase transition for fg = 0.5 is also illustrated by the dotted line. 
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FIG. 9: Isotherms for a phase transition at T = lOMeV and the initial condition /s = 0.15. The 
sohd Hne AC is built up by Maxwell construction, and the system evolves along AC. 
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FIG. 10: Entropy as a function of temperature at constant pressure p = O.lTMeVfm ^ for 
strangeness fraction fg = 0.50. The entropy evolves continuously through the phase transition. 
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FIG. 11: Heat capacity as a function of temperature at fixed pressure p = 0.17 MeVfm~^ for 
strangeness fraction fg = 0.50. The discontinuity makes a strong proof for the second order phase 
transition. Note the logarithmic scale used on the left. 
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